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"fj Abstract. We investigate up to the fourth order normalized factorial moments of 

^) free-propagating and pulsed single photons displaced in phase space in a phase-averaged 

manner. Due to their loss independence, these moments offer expedient methods for 
quantum optical state characterization. We examine quantum features of the prepared 
displaced states, retrieve information on their photon-number content and study the 
reliability of the state reconstruction method used. 
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1. Introduction 

The observation of quantum hght in phase space is typically associated with honiodyne 



■^ detection [I]. However, photon counting also provides attractive approaches for 



cn studying the distinct quantum character of the quasi-probability distributions |2||3]. 

^ Using the photon-parity operator, for example, even individual points of the Wigner 

^ function can be directly probed despite the Heisenberg's uncertainty relation of the field 

• ^ quadratures [4]-[6]. In the optical regime, this technique has already been utilized for 

/\ the characterization of cavity confined optical states (7 10 as well as free propagating 

H — 



d light fields 11 14 . Even though the photon-number parity can be straightforwardly 

deduced from the measured photon statistics, the detection losses of practical photon 
counters tend to destroy the real quantum characteristics featured in the photon-number 
distribution. Therefore, sophisticated methods for testing the nonclassical character of 



the loss-degraded states via their photon-number content have been developed 15 , 16 



and techniques allowing a loss-tolerant reconstruction of photon statistics are of great 



importance 17-22 . Still, an accurate determination of the state's properties from the 



loss-degraded data is challenging, especially when the studied state incorporates higher- 



photon number contributions 23-25 . 



Genuine quantum features can fortunately be recognized even without access to 
the state's complete phase-space representation. One of the pioneering techniques 
was introduced by Hanbury Brown and Twiss 26 , and in an extended form 27 



their experiment allows us to access the higher order factorial moments of photon 
In general the m-th order normalized factorial moment is determined 



number 28 30 

as 
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(m) 



(n(-))/(nr, 

where {n^"^^) = {: ra™ :) = X]n'^(^ "~ ^ 
evaluated either as normally ordered ( 
or via the photon statistics g{n) with n being the photon number 31 



(1) 

. {n — m + 1) g{n) and {n) = {n^^'') can be 
moments of the photon number operator n 

Usually, the 



normalized form of the factorial moments can be extracted loss independently 29 32 



However, care has to be taken when detecting multimode states since the individual 
modes may suffer from different amounts of losses 33 . Nonetheless, these moments 



provide versatile alternatives for investigating the properties of quantum optical states. 

Regarding single photons 34 , already their second order normalized factorial 

moment, which ideally takes the value g^'^^ = 0, can be employed as a valuable 

characterization tool. The real measured values are widely used to classify the practical 

More generally, however, observing g^"^"* < 1 can be regarded 
and as an indication of the sub-Poissonian 



35 



single-photon sources 

as a signature of the nonclassicality 36 



photon-number characteristics 37 39 . Nonetheless, when single photons are displaced 



their g^'^'' values gradually increase and finally exceed unity, which signalizes a super- 
Poissonian photon-number distribution 40 . In this region also a more sophisticated 



test for the nonclassicality should be found 41 , 42 . Further, once having accessed 



the photon-number content of displaced single photons 14 , 43 many other intriguing 



phase-space features can be directly scrutinized such as the nonclassical oscillations 
in the photon statistics 44 , 45 . Moreover, the factorial moments of displaced states 



indeed provide routes for accessing more elaborate moments of the photon creation and 



annihilation operators 27 



Here, we measure up to the fourth order normalized factorial moments of phase- 
averaged and displaced single photons directly in a coincidence counting experiment 



by employing the time- multiplexed detector (TMD) 146, 47 that has proven to be 



a powerful tool for measuring the higher order moments of pulsed quantum states 
of light |29i]. Even without access to the complete photon statistics we can loss- 
independently observe quantum features in the prepared states. Further, by calibrating 
the mean photon number of the loss-degraded states we gain information about their 
photon- number content. At low detection efficiencies, it is generally a highly nontrivial 
task to invert the action of losses 48 . Further, in order to estimate the reliability of 



the state reconstruction often numerical methods are applied such that the effects of 
statistical fluctuations, losses, and highest resolved photon number can also be taken 
into account 14 , 25 . In contrast to the ordinary loss inversion, our technique allows us 



to study the effect of different experimental limitations separately from each other. As 



a consequence, we can directly define the boundaries that an experimental realization 
sets for the state reconstruction in phase space. Further, the artifacts introduced by the 
nonideal detection can be recognized in a straightforward manner. 

This paper is organized as follows. In Sec. |2] we review the properties of factorial 
moments, which are connected to the photon statistics via the moment generating 
function. In Sec. |3] we survey the properties of displaced single photons and review 
the effects caused by experimental imperfections. In Sec. |4] we utilize the factorial 
moments to investigate the characteristics of the prepared displaced single photons. 

2. Retrieving state characteristics via factorial moments 

The normalized factorial moments can be directly utilized for different state 
characterization tasks such as discrimination between sub- {g^"^^ < 1) and super- 



Poissonian {g^'^^ > 1) photon-number distributions 49,50 or classification of different 



quantum states [29]. Additionally, they provide several alternatives for the direct 



examination of quantum features 29,51 . One option for investigating phase- insensitive 
nonclassical behavior in a single mode is to study whether it is possible to violate the 
criterion 
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(m+l) 



> g'^""^ > 1 (2) 



that classical states obey 52 , and we treat the multimode states equivalently. For this 
purpose, also the measurement of the normalized factorial moments with orders higher 
than two becomes relevant when regarding states with super-Poissonian characteristics 
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A more detailed investigation of the state's photon-number content via the factorial 



moments is possible by employing the moment generating function 32 , 53 that is 
described in terms of a real valued variable < /i < 2 as 

M(/i) = ^f,H(l-/i)". (3) 

n 

At the upper bound, that is when /i = 2, the moment generating function provides 
information on the photon-number parity. Nevertheless, equally interesting is the lower 
bound since the derivatives of Eq. ([3]) at /x = are directly connected with the factorial 
moments by 



("""')= -x: A#WL.o- (4) 



Once having accessed the factorial moments the expression in Eq. (|3| can be re- 
written as an expansion 



and the photon statistics is again reconstructed from the moment generating function 
in Eq. ^ by 

^ ^^ m>n ^ ' 

which obeys the conventional normahzation X^n ^('^) ~ -'-• Thus, the photon statistics 
is gained by summing up the factorial moments with proper weight factors. However, 
we note that this method for reconstructing the photon statistics can only be successful 
when the expansion in Eq. ([s]) converges near /x = 1. 

The loss tolerance in Eq. (pj) is achieved after re-writing (n*-*"^) = g'^'^^ ('^)™ 
and deducing the mean photon number (n) from the loss-degraded measurement via 



(^)iossy ~ '^ ('^)' ^^ which T] is the detection efficiency 32 . Thus, apart from recording the 
different orders of g'^"^\ this method further entails a calibration of r/ and measurement 
of (^) lossy- ^^ ^'^^ mean photon number can be calibrated accurately, the reconstruction 
of photon statistics via normalized factorial moments becomes especially expedient at 
low detection efficiencies. The only constraint lies in measuring enough orders of g^^^ 
with good precision during a finite integration time. As a consequence, the highest 
statistically accessible moment, which is known from the experimental data, limits the 
possibilities to completely reconstruct the state's photon-number content. Natural limits 
are determined by the physical bounds < qin) < 1. More stringent conditions may 
be found by investigating the distinctive features in the photon statistics of the studied 
states. 

3. Modeling free-propagating displaced single photons 

Even though displaced single photons have been studied in several experiments 
[8| [l0|[l4| , [43] , the generated states are seldom ideal and imperfections in the preparation 
process have to be taken into account. We first investigate the properties of ideal 
displaced single photons and then regard the effect of experimental imperfections. We 
consider a displacement with a mode mismatch and take into account higher photon- 
number contributions of the prepared state. 

In the single-mode picture, an ideal displaced single-photon state is described as 
D{a) |1), where -D(a) [with the hermitian conjugate Z)^(a)] defines the displacement 
by an amount of a and |1) is the single-photon Fock state. We deduce its factorial 
moments by evaluating the normally ordered mean values in Eq. ([I]) with the help of 
the transformations D''{a)aD{a) = a + a and D'^{a)a'^D{a) = a^ + a* of the photon 
annihilation (a) and creation (d^) operators. After a straightforward calculation, the 
normalized factorial moments of the ideal displaced single-photon state can be expressed 
as 






(m) _ |"| ^ '{III' -r \tx\ ) ,-. 

^ideal — f^ , |„,|2Nm ' y' > 



These moments as depicted in Fig. [IVa) rapidly grow from zero when increasing the mean 
photon number given by {n)-^^^^^ = 1 -|- |ap. By following the behavior of the second 
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Figure 1. Four orders of the normalized factorial moments, shortly 5*^™^ {m = 
2,3,4,5), with respect to the mean photon number (a) for the ideal displaced single- 
photon states evaluated according to Eq. (l7| and (b) in the case of no mode overlap 
between the single-photon Fock state and the reference state as predicted by Eq. (|9|. 
Dashed lines separate the sub- and super-Poissonian regions. 



normalized moment, one directly concludes the gradual transition in the photon-number 



characteristics of the displaced states 40 , 45 , and this moment reaches the maximal 
value of f^itjeai ~ 1.333 at the displacement |a|max = V^- Further, the ideal displaced 
single-photon states always violate the criterion in Eq. ([2]). However, as seen in Fig.jIVa), 
the verification becomes more and more challenging when increasing the displacement 
since moments with higher and higher orders have to be resolved. 

In practical applications, the displacement can be implemented with an asymmetric 
beam splitter, at which the studied state is overlapped with a coherent reference state 



see e.g. 25 and the references therein). The mismatch of the two overlapping modes in 



temporal, spectral or spatial degrees of freedom can be modeled with a simple overlap 
factor A^. In order to evaluate the required mean values in Eq. ([I|, the photon-number 
operator can be replaced with 

n-^T] hes (8) 



V 



D\VMa)a^aD{VMa) + (1 - M)\a\ 



where rj is the total detection efficiency, Ai is the mode overlap, and a is the amount 
of the applied displacement 25 . The effective photon-number operator Ues in Eq. (|8 



is a sum of a displaced photon-number operator and a background term. 

The total detection efficiency cancels out when evaluating the normalized factorial 
moments, and we can write them loss independently in the form 



(m) 



n!!?,:) 



E. 



'-cflf 



{nC^X (n(™-'=)). 



in which (n^"^-*) 



^ = [(1 — A^)|Q;p]™' describes the properties of the background and 



(9) 



= {D\VMa)d^"'a"''iD{VMa)) (10) 

predicts the behavior of the displaced part. If the mode overlap is imperfectly aligned, 
the results of the measurement change drastically. As shown in Fig. [lib), in the case 
of the single-photon Fock state the effective normalized factorial moments cannot take 
values larger than unity when A^ = 0. Therefore, the super-Poissonian region becomes 
a loss-independent indicator that the displacement takes place, in other words that 
A^ 7^ 0. Nevertheless, our model in Eq. ^ provides inaccurate results if only single- 
photon Fock states are regarded. This can be as corrected by considering the higher 
photon-number contributions of the real single-photon source. We assume that the 
single photon is prepared into a photon-number mixed state [54], and we take the 
higher photon-number contributions of it into account in Eq. ( 10 ) when fitting Eq. ^ 



against the values measured for the displaced single-photon states. The overlap factor 
Ai is held as fitting parameter. Moreover, the amount of displacement applied to 
the prepared single photon can be straightforwardly extracted in our model from the 
effective mean photon number. Even in the case of imperfect mode overlap this is given 
by {fies) = (^)sp + |q^P) i^ which (n) is the mean photon number of the prepared single 
photon. 

4. Experimental investigation of displaced single photons 

In our experiment, shown in Fig.|2| we heralded single photons from a pulsed waveguided 
twin-beam source based on parametric downcoversion. For preparing displaced states, 
the heralded single photons were overlapped at an asymmetric beam splitter with 
coherent reference states in a phase-averaged manner. The amount of displacement was 
controlled by altering the mean photon number in the reference beam. The heralded 
displaced single photons were then coupled to a TMD for detection. Our experimental 
arrangement is similar to the one used [l4] except that the TMD is employed for 
measuring manyfold coincidence counts and not the so-called click statistics. We note, 
however, that the values of the normalized factorial moments could also be deduced 
via the loss-degraded photon statistics that can be retrieved from the click statistics 
recorded with TMD by taking into account the intrinsic detector characteristics I50|. 



As described in 29 , the higher order correlations in a light beam can be accessed 
with the TMD detection scheme, and the time-integrated measurement delivers the 
desired expectation values J55J. Thence, the value of g^™"^ is extracted from the raw 
data by dividing the probability to measure a coincidence click between m selected 
temporal TMD bins by the product of the single click probabilities in these bins — all 
of them conditioned on the detection of the herald. In order to apply the measured 
normalized moments for the loss-tolerant reconstruction of photon statistics, we further 
require a calibration of the detection efficiency such that the mean photon number 
can be determined. The perfect photon-number correlation between the twin beams 



empowers us to estimate the detection efficiency according to Klyshko 56 . For this 



purpose we block the reference beam, subtract the amount of accidental counts from 
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Figure 2. Experimental setup in accordance with [M]. Ti:Sapphire laser pulses 
(wavelength 796 nm, bandwidth 10 nm, repetition rate 2.4 MHz) are frequency doubled 
in a nonlinear crystal (NC) and coupled to a 1.45 mm long, periodically poled, type-II 
KTiOP04 waveguide (WG). While color glass filters (GF) block the residual beams, 
in the WG generated twin beams are separated at a polarizing beam splitter (PBS). 
Selected by the orientation of a half- wave plate (HWP), one of them — the herald — is 
filtered to a bandwidth of 1 nm with an interference filter (IF) and sent to a silicon 
avalanche photodiode (Si-APD). Meanwhile the heralded single photon is coupled to 
the same spatial and temporal mode with a cross-polarized coherent reference beam 
attenuated to the single-photon level with a neutral density filter (ND) and filtered 
to a bandwidth of 1 nm. The cross-polarized beams are then sent through a spectral 
filter with a bandwidth of 1 nm. The displacement is realized in a phase-averaged 
manner with a PBS placed after a HWP, whose axis was minimally tilted from the 
polarization direction of the heralded state. Finally, the TMD divides the incoming 
light pulse in three subsequent symmetric beam splitters (BS) into eight temporal bins 
that were detected with two Si-APDs. Components in dashed boxes were implemented 
with fiber-integrated optics. 



the coincidences between the twin beam detection, and compare this number to the 
amount of single counts in herald. The mean photon number of the displaced state 
merely follows from dividing the first ttnnormalized factorial moment — the probability 
of measuring a single click conditioned on the detection of the herald — by the estimated 
detection efficiency. 

We first study the values of the moments g^'^^ to g^'^^ in order to categorize the 
photon-number properties of the prepared states. For the prepared heralded single- 
photon state we extracted the values g^'^^ = 0.184(4) and g^^^ = 0.04(2), the accuracy 
of which is limited by the statistical fiuctuations only. In contrast to a genuine 
single-photon state, one clearly recognizes an additional two-photon contribution in the 
heralded state. This is the trade-off from a rather high pump power — on average 12 |J.W, 
which nevertheless yielded a heralding rate of 11.8 kHz. When the prepared single 
photon is now displaced, the values of g^'^^ to g^^^ gradually increase. This behavior is 
shown in Figs|3](a-b) with respect to the calibrated mean photon number. Maximally we 
observed the values f^inax = 1.148(6), gralx = 1.18(3) and g^nlx = 1.00(15). A comparison 
of these results with the ones gained for the case of vanishing mode overlap [Figs. IsFc- 
d)] reveals, as expected, the rapid appearance of the higher order normalized factorial 
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Figure 3. The measured values (symbols) of (a) g*^^' and 5*^'^-' as well as (b) g^'^^ 
in dependence of the calibrated mean photon number when the mode overlap was 
optimized to A^ = 0.71(2). (c) The same as panel (a); and (d) the same as panel (b) 
for the case of deliberately mismatched mode overlap. Solid lines are fitted, whereas 
dotted lines predict the results for A^ = 0. Vertical errorbars are given by statistical 
fluctuations only, whereas the horizontal ones are dominated by the accuracy of the 
efficiency estimation. Dashed lines are a guide to the eyes. 



moments, when the prepared single photon is displaced. In contrast to our earlier 
studies 14 , we now fit our model in Eq. (|9| against the properties of the higher photon- 
number contributions. Our model assumes that all the photon-number components of 
the prepared single photon possess the same overlap factor and a reasonable match 
for the data in Figs |3](a-b) is found. If the overlap factor was drastically lower for 
the higher photon-number contributions — in our case the two-photon contribution — 
we would expect to encounter difficulties when fitting our model against the values 
measured for the third or higher order normalized factorial moments. 

Even without information on the complete photon statistics, we conclude from our 
results in Fig. |3]^a) that the expected transition between the sub- and super-Poissonian 
photon statistics takes place. Further, our displaced states show quantum features at 
small values of displacement. As seen in Fig. p[a), the inequality g^'^^ ^ 1 guarantees the 
nonclassicality of the prepared states below the mean photon number of approximately 
1.6. Our results in Fig. k3[a) clearly also violate the classicality via g^^^ ^ g^'^\ However, 
this criterion becomes inadequate close to the mean photon number of 1.9, beyond which 
we are unable to detect quantum features in the prepared displaced states. In order to 
do so, a more accurate measurement of the fourth normalized moment g^^^ than the one 
shown in Fig. [stb) is required. 

Next, we examine, how much information on the displaced state's photon-number 
content can be deduced via the measured moments. We estimate the photon-number 
contributions of the heralded and displaced states by plugging the measured values in 
Figs |3](a-b) into Eq. ^ together with the calibrated mean photon numbers. For the 
heralded single photon, the photon statistics of which is shown Fig. |4Fa), we extracted 
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Figure 4. The reconstructed photon statistics of the (a) heralded and (b) minimally 
displaced single-photon states. 



the loss-calibrated mean photon number of (^)sp = 1.07(3) that was estimated in a 
measurement with little less than one percent detection efficiency per TMD bin. Further, 
we note that the accuracy at which the mean photon number is deduced depends not 
only on the statistical fluctuations but also on the approximative estimation of the 
detection efficiency. As shown in Fig. |4Fb), we can recover the behavior expected for 
the displaced single photon — increasing vacuum contribution and decreasing one-photon 
contribution — at the mean photon number of 1.30(4). 

However, the limited resolution restricts the acceptable phase-space displacement 
and sets a boundary to a region, in which the photon statistics of the prepared displaced 
state can reliably be reconstructed. Being able to statistically resolve moments only 
up to the fourth order, in other words the experiment delivers (y((">^) = 0, the two 
highest accessible photon-number contributions are deduced according to Eq. ^ as 
^(4) = (n)V4! g^^^ and ^(3) = {nf /{3)\ [g^^^ - g^^^n)]. The non-negativity of 
the photon-number components provides us with a condition (n) < g^^^g^^^ for the 
accessible reconstruction range. By employing the fits in Figs |3](a-b) we get the bound 
(n) < 1.4. Looking closer at the individual photon-number contributions in Fig. M 
a second, stricter condition is obtained by studying the boundary, close to which 
the reconstructed photon-number contributions start to significantly deviate from the 
expected behavior. As the highest resolvable photon-number component, in our case 
^(4) [Fig. [5]^e)-dashed hne] , gradually increases with respect to mean photon number, 
it eventually surpasses the contribution g{3) [Fig. |5](d)-dashed line]. This is an artifact 
not expected in the photon statistics of the displaced single photon at the studied 
region. Therefore, we use the limitation ^(3) > ^(4) as the boundary of the reliable 
reconstruction range for the displaced prepared states and gain the ultimatum (n) < 1.3. 
Furthermore, we note that the vacuum component [Fig. |5](a)-symbols] of the prepared 
displaced states could be fairly reliably extracted from our measurement even beyond 
the determined range, whereas the higher photon-number contributions [Figs [stb-e)- 
symbols] are not as resilient. 

In order to accurately reconstruct the complete photon statistics outside the 
determined region, one is required to resolve moments with orders higher than four. 
Thus, we predict the values of the moments g^^^ and g^^^ in order to interpret their 
influence on the photon statistics. As depicted in Fig. [5] the accessible reconstruction 
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Figure 5. The reconstructed photon number contributions e(0)-£i(4) (a-e) for the 
displaced states prepared in Figs Isla-b). Symbols originate from measured values. 
The (orange) dashed lines present the truncation of moments up to g^^^^ and is gained 
from the fits in FigslsFa-b), whereas the (black) dash-dotted and dotted lines include 
the prediction of moments up to g^^' and g^^' , respectively. The (gray) solid line 
represents the theoretically expected photon statistics. 



range gradually increases when moments with higher orders are resolved. Clearly, the 
method used provides us with insight into the reliability of the state characterization 
at low detection efficiencies. In our example it allows us to estimate the bounds 
of the reconstruction range in phase space and to recognize the artifacts caused by 
the nonideal detection. Systematical deviations in the reconstructed photon statistics 
become apparent if enough orders of the normalized factorial moments cannot be 
measured due to the limited detection time. Further, the statistical fluctuations in the 
measured moments and the precision at which the mean photon number is extracted 
affect the accuracy of the reconstructed photon statistics. Fortunately, their effect can 
be investigated separately from the systematical deviations and in our case the errorbars 
are dominated by the fluctuations in the resolved moments rather than by the inaccuracy 
in the mean photon number. In summary, our results can give direct specifications for 
the experimental parameters when studying the fine structure in the photon statistics 
of displaced single photons at low detection efficiencies. 



5. Conclusions 

We measured up to the fourth order normalized factorial moments of displaced single 
photons in a loss-independent manner. By studying the second normalized moment, we 
con&med as expected that the sub-Poissonian photon-number distribution of a single 
photon gradually moves toward the super-Poissonian photon statistics when the state 
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is displaced. The prepared displaced states further violated the classicality even in the 
super-Poissonian regime. Moreover, the measured moments provide means for the loss- 
tolerant reconstruction of photon statistics after determining the mean photon number 
of the studied state. However, it is essential to accurately measure enough orders 
of the factorial moments in order to reliably reconstruct the investigated properties. 
Our results show the versatility of the factorial moments for the state characterization, 
and we believe they prove to be useful in examining genuine quantum features at low 
detection efficiencies. 
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